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ABSTRACT 

The radiation impedance of a finite array of finite cylinders is 
determined by a boundary value problem technique. Tho final equations 
for the radiation impedance are stated without rostriotion on cylinder 
length or radius, oylindor spacing, or on the amplitude and phase of 
the radial velocity of the cylinders. The derivation is made for the 
case of each cylinder in tho array having a uniform radial velocity over 
a sector of the cylinder, but tho method of solution is not limited to 
such a velocity distriLutione The results of a numorical calculation of 
the effeot of eloctrical beam steoring on the radiation impedance of a 
specific array are shown to illustrate such effect and the uso of the 
derived equation for radiation impedance. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by lire John Hickman of the U. Se Navy Llectronics 
Laboratory and Professor Iie Medwin of the U. S. Naval Postgraduate School 


in this investigation. 
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= Wave velocity 
e Base of natural logarithms 
j \ 1 
Ae aes 
k Wavelength constant, k = 
p Instantaneous acoustic pressure 
ne Wavelength 
a) Volume density 
CD 


Angular frequency 


The notation associated with the Bessel Functions used in this paper 
is taken from N. W. McLachlan, Bessel Functions for Engineers, Oxford 


at the Clarendon Press, 1955, 


J,(z) Bessel function of the first kind 
bez) Bessel funotion of the second kind, as defined by Weber 
K (2) Modified Bessel function of the second kind 
(1), (2) | 
Ay Bessel function of the third kind. 
HAs) = J 
. tz) (2) # jY (2) 
(2) 
He mee, (2 ) 


— 1 
\ The positive square root of a real, positive quantity. 





1. Introduction 

The radiation impedance of a finite array of finite cylinders in an 
infinitely long, rigid, sylindrical baffle has been investigated by D. He 
Robeye> Robey considered the problem of uniformly pulsating cylinders 
and his final conclusions are based on a cylindrical array of radius 
negligible in comparison to the wavelength of the radiated sound. It is 
the intent of this paper to investigate the radiation impedance of such a 
finite array of finite cylinders with fewer restrictions imposed on the 
radial velocity distribution and size of the array. It is also the intent 
of this paper to apply the results of such an investigation to the study 
of the effects of electrical beam steering on the radiation impedance of 
a finite array of finite cylinders. 

The radial velocity distribution on each pulsating cylinder will be 
specified as a boundary condition to the wave equation which is to be 
solved. However this is not enough to completely specify a unique solu- 
tion to the wave equation. Therefore it is proposed that as the radial 
distance from the cylindrical radiator increases, vanishing boundary 
values for pressure and velocity aro required. This condition again 
proves too weak to determine a unigue solution. Therefore a more stringent 
boundary condition in addition to the condition of vanishing boundary 
values is required, ie. "The Sommerfeld Condition of Radiation.” That is, 


"The energy which is radiated from the sources must scatter to infinity; 


De H. Robey, On the Radiation Impedance of an Array of Finite Cylinders, 


Je Acouste Soo. Ame,27, pp706-711, 1955 


2A, Sommerfeld, Partial Differential Equations in Physics, Academic 
Press Ince, 1949, p 189 





no energy may be radiated from infinity into the prescribed singularities 
of the field (plane waves are excluded since for them even the condition 
ue O fails to hold at infinity)." 


The wave equation which is to be satisfied is stated in cylindrical 


coordinates: 
27g) a ody a 2b y 2b >t 2g 
aa ig ge) Se r” Bias ee ee a+” 

The scalar velocity potential function is understood to be £(r,e,z,t). 

The z axis of the cylindrical system corresponds to the longitudinal axis 

of the cylindrical array (see figuro 1). 

The driving forces which are to be imposed on the elements of the 
cylindrical array will vary harmonically with time. Therefore for the 
results of this paper to conform to standard electrical enginecring 
practice of assigning positive reactance to inductanco, a time dependence 


of the form eJ“% will be assumed. That is Z(r, e,2z,t) = B(r,e,z) ed#*, 


After substitution for B(r,ez,t), the wave equation is reduced to the 


following: 
So yk 42d et 2 yes: ee a 
r OZ c 


A formal solution to the above equation is discussed in Appendix l. 
The radiation impedance of an array of finite cylinders will now be 


determined. 





20 Derivation 

The first array to be investigated is the least restricted in 
boundary conditions of the arrays to be considered. The array is to be 
constructed as follows (see figure 1): 

1. There are m active cylinders of length = 2f., radius = a, 


and center to center spacing = 2b. 


2e The cylinders are mounted coaxially in an infinitely lone, 
rigid, cylindrical baffle of radius = a. 

The pulsation of any active cylinder is to be described as follows: 
le On the given cylinder, the radial velocity is constant 
over the length of the cylinder for a given aximuth angle. 
2e On all cylinders, the radial velocity is uniform between azimuth 
coordinates -& to #a and zero elsewhere. 

Therefore the radial velocity of the jth cylinder may be described by the 


following Fourier Series: 


Oo 
Me get 
Vv, Ve eJ AL cos ne (2501) 
n= Q 
Where A, = &, Sin ne » yn al, ne O 
no 


en = ae is Near cre 


and Wi = the phase of the pulsations. 
The scalar velocity potential Z(r,e,z,t) must satisfy the following 


wave equations 


2 
Vi gs = 2 p (2.02) 





The scalar velocity potential is assumed to vanish as r-» oo and the 
Sommerfeld Radiation Condition is assumed to be valid with respect to the 
radial coordinate. 

The time dependence of the scalar velocity potential is taken to be 
O(r,0,z,t) = P(r,e,z)eJ” as discussed in Section I, The scalar velocity 
potential O(r,e,z) must then satisfy the Helmholtz Wlavo Equation 
(Vedx%) B= 0. (2.02) An eigenfunction which is a solution to the 
Helmholtz Wave #quation, satisfies the vanishing boundary conditions and 
the Sommerfeld Radiation Condition, is developed in Appendix l. 


This eigenfunction is stated: 
sin(ne) , =j27 fz (2) 


Bry (r,e,2) = ae cos(ne) Ba VA) (2.03) 
Where fe a k* (Bae (2.04) 
fez x2 - (27 7 )@ (2 rae neato 
Aa -j\ary)? = x? (2p) > x? 


It is now proposed to determine a @ (r,e,z), constructed from a sum 
of the eigenfunctions Pn, ¢ (r,e,z), which will satisfy the boundary condi- 
tions at the surface of the cylindrical array. The technique to be 
employed in determining the complete scalar velocity function is a 
logical development from the technique of L. Ve King® used in the solution 
of the radiation impedance of a piston in an infinite rigid baffle. The 
constant Any transforms into a fe the eigenvalue f becomes a 
continuous variable. 


The sum of eigenfunctions is taken to be: 
OO 


~ (2) 27 
B(r,e,z) oe el Cs) cos(ne) H, (Ar) elev Py z ay 


n20 \ 
U 
-~OO0 


205.) 


ere. Kine weommene Acoustic Radiation of the Pilex6-Electric 
Oscillator and the Effect of Viscosity on Transmission, Can wJourn, Res., ll, 
pp 155=155, July-Dec.1934 


° 4. 





The velocity in the radial direction is: 


vy(r,e2z) = - Se g (r,e,z) 


00 
= (2) (2) 
vi(r,a,2) \ ei(p) Saath er) near) cos(ne) od@ be 
n=0 £ 


-00 


(2.06) 


The radial velocity distribution of the finite array will now be 
investigated. From the boundary conditions as defined, it is possible to 


write an expression for the radial velocity amplitude at the surface of the 


array: 


v,(a,0,2) = 


ore) 
a) cos(ne) 
=0 


n 


(2.07) 


The function g(z) is equal to zero on the rigid baffle and on the ith 


cylinder is equal to Vj eJ?i, If the two expressions for the radial 
velocity, (2.06) and (2.07), are compared; the equality of the two 


quantities is possible if they are equal term by term, that is: 


: 7 . (2) (2) ie 
g (2 )An es \st59 Ee eje” Pf Z ay 


(2.08 ) 
Define a new quantity, (2) (2) 
G(f) = “5 £.(F) ok Ufa) = Hema (a) (2.09) 


i) 





Then 


00 
g(z) = G(y)e Ai2r¥2 ay (2.10) 
=o 
It can now be seen that G(f ) and g(z) are a Fourior Transform Pair. 


Therefore: 


G(Y) = e(z) e"J27F 2 ay (2611) 


-00 


bef b¢A 
: \ (0) ew7J@7P 2a, | v, od¥1 oo I7™F 2a, 


~0o bef 
(2i-1)b~HZ 
wh ae weer ao vi eo ie J27P 2 ap ¥ 
(2i-1 )b-2 
(2m=-1)b4#A re) 
Lips earn ae Vneje me “i27P 2 ae o (0) e727 F2 a, 
(2m-1 )b=A (2-1 be 


sin 272 
a mA 


Pues eae es (21-1 Je ane D, 


G( YP) = V; sate 1,-jer F b 


f covet Vin ol hm e-J27¥ (2n-1)b — 


G(f)s y, oi V4 grje27¥ (2i-1)b , sin a 
: vr 
Fe (263281 
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It is now established, after substitution: 


™m 
2A, sin(2mpB ) > Vi Bays e7Je ry (2i-1)b 


i=] 


2 L 
perf [Badaa) = Har Ha)] 





(2.10) 
and 
[ 90 
~ 2A, cos(ne) sin(27 pe yz, (?? r) 
f(r,e,z)= >! Le : if - x 
n=O} Tia tea fea) _ ae K0)| 
[-oo" 
7 


m 


= S Vied¥i 9 ~j27F (2i-l)b ej27Y z dy 
Gol 
Ts 
eee 


To obtain the average pressure, one uses the force equation 


Ve 4@q 70 
Since V = - VG and for the condition of time dependence eJ¥t | the 
force equation takes the form: 
p(r,e,z) = jv ofA (r,e,z) 
The averazte pressure amplitude on the active area of the nea cylinder is 
defined as By. 


+o (2h-1 )b+p 
Phe ror B(a,e,z) de dz 


-% (2h=1 )b=p (2.15) 





The averaging process is oarried out and the expression for the 


average pressure is, upon substitution for A,: 


ona) 
"2 
GO 9 ”, H ( wa) 
a: \ sin“ (nx) penal ¢ 
Przjro / A ° (naj (ame? % a ) - H, (sa) 
Y ajJ- & 
n20 P pe Pada ¥ n=l, } 
-00 
Oo 
- > gihie 34 7 F- (i-h)b ay 
i=] . 
(2.16) 
The average radiation impedance per unit area of the active area 
of tho h&® oylinder will be defined as <z> n° . 
(2), = bh 7 
Vy ett 
(eli 
i 
Further def ae 
6r eiines G" a ‘ Po 
itercae 2 (Ay ae 
= j- Wap 2-2 (ay Boa 
(2.18) 
: \ 
Then the radiation impedance per unit area of the active area of the 
pa cylinder may be written: 
ore 
2 ae (2) 
" % B sin (na) sin” (2 7r Jeeay, 0" ka) z 
ead ee a ae 2 ( | 
a (n) (27 ¢ Pp )*r [tay (a xa) 22) ra) 
n= 
-00 





m 
| 
. <a 5 Yi -Vn) Sola R (i-h)b ay | 

= Vh | 

i-0 
= (2alep 

The radiation impedance of an array of finite cylinders which pulsate 
uniformly with arbitrary amplitude and phase over a sector of each cylin-~ 
der has been determined. Zvaluation of equation (2.19) involves numerical 
intecration of an infinite integral. The integrand of (2.19) must be 
carefully examined in the region of i> =. A further handicap to 
evaluation of (2.19) is the fact that second order Hankel functions of 
complex arguments must be evaluated. These qualities of (2.19) are 
especially undesirable if tho radiation impedance is to be evaluated by 
digital computer techniques. Equation (2.19) will now be reduced to a 
form suitable for numerical integration. 
The following relation is stated and used in deriving a new form 
of equation (2.19).4 
Ky(jx) # = de d5™ (vol) 4 (2) (x) 
1 


<8 


e Sha (2) 
sin® (na) sin“ (na) sin’ (27 A (gvka ) _ 


(na) 


n=O (2 ir PB )Eeg [ees (ea Joh, 3 (ee) | 


on re 


Vij (Vi-~h) . 
<> - cos [ 4% Y b(ieh) | ay + 
i=l 


—_ 


aN McLachlan, Bessel Functions for Engineers, Oxford at the Vlarendon 
Press, 1955, ps 204, 





sin"(2% ye )(-j) Kn (jo ka) - 


(2% $72)? © fxr (joka) pK, 7 (Joke) 
; 


exe) 
Oo 
> en sin“(n) 
(n ) 
Lv 
oe 
m 
ve Saree 
a oi Hi pe 30s | 47 p b(i-h) | slag 
Vn -_ 
i=] 
(2.20) 


The integrals will now be considered separately. Let 


(2) 
- sin eS 308 [47 ¥ (i-h)b | dy 





= 
(2% FB )e ¢ hoy ( ~ka) - a (ea) 
(nen) 
Let APs sin § 
then when f= 0 ; Fe 0 
i 25 ; f | +7 
Cara a ware Ne = cos § 
dvs cOs ¢ a f 
a 
Therefore: 
Lr 
sin® sin In (*) & GOS 3 
a KP — Paes) = cos | 2kb(i-h)sinf]d § 
F (kf sin¢ ) Hn A) (ka cos ~—) = 12) (ice cos f ) e 
(2eo2F 


10 





Consider: 


ore) 
sin (eee 3) (<4) Knlig ka 
I, = _ aes COs [airy (i-n)b] d¥ 
(Qumae eg [Knoft (a ka ) A Kyi (jo Ka)] 
aL J 
a 
teece) 
Let Af s ose § 
Then when a ; f = arr 
x 
= 00 ; f = 0 
) a | 
2 -j yar) -1 = -j cotf 
qdfve - osc § cot: — 
Therefore: 
ue 
oe K_ (ka cot — ) 
Ie \| = i pieeaa) eee) ee cos [ 2kb(i-h Josof] as 
(kB )* csof — Kyfy (ea cot £) f Kyy (ka cots ) az 
O 
(2324) 


The expressions (2.22) and (2.24) are substituted in (2.20) and the 
resulting equation for radiation impedance per unit area z jj, takes the 


more tractable form: 


ove) (2) 
# >| v2 H (ka cos ?) 
Aaa i800 <A sin’ a a sin (kp sin ¢) n 5) = 


(now )* (kp sin §) HA Ua gost )=Hn 4 (lea cosé) 


Va 





m 
ee 
>) peer? cos [2ieb(i-h )sinf] a § a 


‘cies, | 

: az 
race) 8 9, SM, i sin® pe so —Kn(ka cot f) x 
(nx) : (xB) csc 8 Ky fy (ka cot t) A Ky2j (ira cot} ) 





vA ers -¢n) cos [2kb (ih oso ¢ | df 


(eco) 

The evaluation of ¢z2)} has been reduced to that of evaluating finite 
integrals and functions of positive real variaovles for which mathematical 
tables are readily available. A further manipulation of (2.25) will be 
made in order to separate the real and imaginary components of <z)> pe 

Let 


1 (2) (xe COS 
Bn (ka cos fF ) en (ka cos § je i (k a) 


2) 2) ‘ 
Ty dilka cos f ) - Hy.y(ka cos¢ ) 
(2.26) 


Then the expression for <2) }, transforms into the following expression: 


12 





fT 


(i2e°c) 











( } 409 6x) ¥ 


wale ay) Tm y A ( $ 400 2) Th Yy }oso (ft) ( ua) =e v 

ee. : T oes fj 2 a df 

3p | fo (4 )axg| 00 (fo *f) 00 ae (Neu) urs i < ox 
w re 


é 


0 
; Ts} f 7 
a Us (< Urs JX) (ou) 
3p | Sure(uesyarg] eoa(d-'4) - ( 2 #00 1)"p. | 200 F, < } 200 ex) € - } Coaygare = : 
= oo 











dp | Jos0(y-t)qx2| soo(4-! Auts - a. (3300 wt) Sn A( - a a sere = \ ‘ie M9 re ol 
5 | ($300 wx) O80 FH) urs (ou) ute res 
w ag 00 
0 


Trt ( Surs 0) 


u 
s(y- 3 -T4) ~ ($200 vx) Qluts 2 sO i ae 
2p uy (y F) Az | ‘oFe) uy A) (3 aR t aA ; 4) BX ) q Cy urs 5) urs (30 &) pars 3 4 =o 


0% CH 


od 9 vor 


oag = 





The radiation impedance per unit area of an array of m finite 
ey lindeven ee pulsate uniformly over a given sector of each cylinder 
with phase and amplitude arbitrary but constant for a given cylinder is 
expressed by equation (2.27). There are no restrictions stated on the 
wavelength of the radiated sound, on the cylinder length or radius, or 
on the cylinder spacing. Although equation (2.27) is for a system with 
uniform velocity and phase across each radiating element, the method 
of solution is not limited to such a velocity distribution. Any 
velocity distribution in azimuth which can be expressed in a Fourier 
Series is possible of solution using such technique. Further any 
velocity distribution in the z coordinate for which a Fourier Transform 
exists is also capable of yiolding a solution. The distribution in the 
z coordinate need not be tho same on each cylinder for a solution could 


be constructed piece by piece for such a system. 


14 





Se Special Cases 

The validity of equation (2.27) will be considered by investigating 
two special cases of the finite array for which solutions already exist. 
The first array to be considered is the single uniformly pulsating 
cylinder in an infinitely long, rigid baffle, Let the array be con- 
structed as follows: 

1. There is one active cylinder of length = 2f and radius = a, 

ee ine cylinder is mounted coaxially in an infinitoly long, rigid, 

cylindrical baffle of radius = a. 
3e The cylinder pulsates uniformly with radial velocity = 
NG Rag 

The Fourier Series equivalent of the third condition as stated above 
is expressable in one terme That is, on the active cylinder the radial 
velocity may be expressed as one term of the Fourier Series (2.01). 


Vr 3 vG Baie 


= 


With reference to equation (2.27), tha above statement or condition 
reduces the summation with respect to n to one term, namely n = O. The 
condition of one cylinder reduces the summation with respect to i to one 
term, namely i= O. Let the radiation impedance per unit area of the 


uniformly pulsating cylinder be <z>. Then from equation (2.27), 
os 


2 


- (eles: / sin” Org sins) Bo (ka cos f ) sin aon ka cos ) dé 4 
ps (kB sin  )° 


raw) 


als 


e 

iS sin” sin“ (kg sin § Z 

7 Eyes g\ at ing) ——~— B,(ka cos é) cos X% (ka cos s) d 
; a. : 


LS 





c 
Wags CSC ) Kolka cot +) : ‘ 
# j8 Be fe ce 
nr (kB ) asc £ 26 ka cot g) 
0 
(50m 
Equation (3.01) is written with substitutions according to (2.26): 
aie 
vz 


us 


Wry. 


en) & 4ae & sin’ (k sin‘ ) Hy '? ica cos & ) 4 


(kp sin § )’ He COs § ) 


T 
ce 
a oi af sin’(k 3cscf) Ko (ka cot é) at 
a=(o A (kp fase § Ky (ka cot : k ) 


io.02) 

Equation (3.02) compares with Robey's equation (18)° for the 
radiation impedance of a single uniformly pulsating cylinder in an 
infinitely long, rigid oylindrical baffle if the variables of integration 
are changed and it is allowed that Robey's equation contains an error in 
signe (This is demonstrated in Appendix II). 

It will now be shown by uso of the Dirac Delta Function that 
equation (3.02) will transform into the well known expression for the 
radiation impedance of an infinitely long, uniformly pulsatins cylinder. 
Equation (3.02) is written in the form of equation (2.19) by trans- 


formation of variables: 


Op, i. Robey, On the Radiation Impedance of Finite Cylinders, J. Acoust. 
Soc. Amo, 27, ppe 706w711, 1955 
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ster e) Ho $2) (~ ka) 5 


<2> = Jes 
; Wei epee, 
- oe a CH, (ka) 
; (sas) 
Consider the expression q 
eee 
lim — sin” (278k) 
217/800 27778 Ce )é 
(3.04) 


Expression (3.04) has the following qualities: 
Was ZOrOr1 Or all ie except Fr = O-« 
2. At JF = 0, (3.04) tends to infinity. 
3e The integral of (3.04) from f equals minus infinity to 
yf equals plus infinity is equal to unity. 


The above is the definition of a Dirac Delta Function. Therefore: 


ee | 
Lim 1, sinterey) - o( P) 





amar Tora CE) 
(5.05) 
and 
OO 
H (2) (a ka) 
lim<2> = jee 2(¥) 4 
2irB-> 00 pe Hy Gis 
-00 
3 (2) (xa) 
lim <2> = jpec 
oN Ra ie Hy (ka) 
(o.070 
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The equation (3.07) which is the radiation impedance per unit aroa 
of an infinitely long, uniformly pulsating cylinder compares with the 
well lmovm expression for such radiation impedance. ®© TRAis concludes the 


discussion of the validity of equation (2.27). 


i Beranek, Acoustic Measurements, John Wiley and Sons, Inc., 1949, 
po 60 
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4. Beam Steering ; 


Tho application of equation (2.27) toéspecific array will now be 
considered in order to study the effect of electrical boam steering on the 


radiation impedance of the arraye The array is to be constructed as follows: 


1. There are nine active cylindors with kf = in, ka = and 
kb = 0.677 
Ze tlhe cylinders are mounted coaxially in an infinitely long, 
rigid cylindrical baffle with ka =7 
3. The cylinders pulsate uniformly over thcir surfaco (symmetry 
in azimuth). 
It is proposed to investigate the effect of electrical beam steering on 
the radiation impedance of the individual cylinders. H#lectrical beam 
steering will be accomplished by phase leading each cylinder a uniform 
number of degrees ( ip) with respect to its position in the array. 


Equation (2.27) is written to conform to the above boundary 


conditions o 
ft 
ie 
<2 > . ean seeaae : iH (2) 
; ieee at sin aeae meee) ee | 1.2 7 (i-h) sing | o .~—sMrcos§) 
of BT sin : 
| tap Se m)roos ) 
ay 
é 9 
saya SS ae : 
nae sin® (s{fc sc f ) od (i-h )y | Th cot E ) 4 
EW )osc § = a [2627 (ichosof] Foor ) 


en) 


IBS 





The result of the evaluation of (4.01) by numerical integration 
using the trapezoidal rule and At = £ is plotted in figures 2 and 3 
for cylindérseaumber one, three, and five. 

ror such an array, the electrical beam steering has but a moderate 
effect on the radiation resistance of any element. [Electrical beam 
steering affects the radiation resistance of the end elements of the array 
to a greater degree than it affects that of the center elements, 

Klectrical beam steering has a more drastic effect on the radiation 
reactance of such an arraye Radiation reactance in this array varies as 


much as ten to one for the center element of the array. 
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APPENDIX 1 


The Helmholtz Wave Squation is restated in cylindrical coordinates: 


a2 1 oo 1 62 a2 ee 
See Ft Set Fro — 


1G 


The solution to the above equation is assumed to exist in a variables- 


separable form: : 
O(r,e,z) = R(r)@(e)Z(z) (A1.02) 


(A102) is substituted into (Al.01) and the resulting equation is: 


A Os) 








jemaR  ,) dk yd ae 2 a*z a os 
R d@r® ° Rr dr @Y® de*°~ Z~  dzé 


o 
It can now be seen that a — is equal to a quantity which is independent 
Z dz 


of z3; therefore the following is true: 


2 
1. — 5 (20 Op (Al.04) 


Where iF is an eigenvalue to be determined by the boundary conditions 
assigned to the array. The solution for Z(z) may now be written. 
Z{z) = By eJeTy¥ z tby qe “jews z (A1.05) 


a: 1 da’ . 
Continuing with the solution, it can be seen that —= is equal 
: @ dee - 


to a quantity that is independont of e; therefore: 


ae os 
4 F = - x! (A1.06) 


Where n is an eigenvalue to be determined by the boundary conditions 
assigned to the array. It will be found more convenient to write the 


solution for 6(e) as follows: 


6(e) = cycos(ne) # d, sin(ne) (Al. 077 
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The assignment of arbitrary eigenvalues to the solutions for 2(z) and @(e) 


determine the form of solution which Rn yp (r) must assume. 

aay, 2 Bay, { [x? _ (278 )? | =} Ray — 0 (A1.08) 
2 ) i = 

ehe r ele r 


In the following discussion, the quantity k* - (2@ )® will be defined as 
equal ton 4, The quantity must reduce to k as (2a7¥ )¢ approaches zero; 
therefore: 

ye fee= (ary? (amy? < 2 

JES aj ((2ar P )2 2 (2m¥)2 7 ke (A1.09) 
To satisfy the Sommerfeld Condition of Radiation, the solution for Ry,» (r) 
must be: 

Rn,p (r) = Ht?) (er) (A110) 

The expression for Rn yy (r) as r approaches infinity is: 
2 gee “p22 22) (Al.11) 


i 
This expression coupled with the time dependence e Jot represents an 





Rn Ge = 


outgoing wavee 

The wave equation has been solved consistent with the Sommerfeld 
Condition of Radiation. It is now possible to state a eigenfunction which 
is a solution to the wave equation; and if the boundary conditions imposed 
on the array are consistent, a sum of the eigenfunctions should satisfy 
the additional boundary conditionse 


The eigenfunction is stated: 


- 


Pay (r,0,2) = An,» 5 ne eee ec ty rr) (Al.12) 


cos (ne 
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APPENDIX II 
To substantiate the statement that equation (18) of Robey's paper is in- 
correct, the development of equation (18) from equation (16) will be inves- 
tigated as well as the validity of equation (16). Further it will be shown 
that equations (19) and (20) of Robey's paper are correct. 
poo 
Ho 62) (Ka ) sin®(f L) 


eS AS <) (Ka) 
0 


Z2qq = j8a pw ay (16) 





2 3, 
Where K = (k°-f 2)& 


Rewrite equation (16) as follows: 


ae (2) 
< = =j fa sin“ (f? L) _ Hitt (Ka) af (A2.01) 


Zaq? eine py K Hy (2) (Ka) 


Where <2qq? = the averase radiation impedance per unit area 
Consider the following expression: 
Pay eins lek (2202 ) 
L>oo TL )* 
The following facts are true: 
1. (A2.02) approaches infinity as f approaches 0. 
2e (A2.02) is zero for all other ¥. 
3. The integral of (A2.02) from f = -00 to FY = “oois equal to unity. 
Mierefore (A2.02) is a Dirac Delta Function, and (A2.01) may be written 


as follows: 
foo 
lim ¢z4,> 2 ~ 30” o(f) 


L-* 66 


H, 6?) (Ka) 
K 11, (2) (Ka) 





ap (A203) 


-OO 
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Thus the radiation impedance per unit area of the uniformly 
pulsating cylinder as derived by Robey reduces to the following as the 


length of the cylinder approaches infinity: 


Hy) (ka ) 


(A2.04) 
un, {?? 


lim)  €Zqq> = ~Jfo° 


L~*> oo (ka ) 


This does not agree with the well known result for an infinite 


(2) 
cylinder. [Further the angle associated with foe) is always in the 


iB] (ka) 
fourth quadrant; thus the radiation resistance givon by (A204) would 


always be negative. 

It oan be assumed that Robey's equation (16) is correot except for 
sign, and it will be shovm that Robey's equation (18) has an error in 
sign but that equations (19) and (20) are correct. Equation (16) is 


restated with the oorrect sign: 


| sin® Cb) Hy”? (Ka) 
Zqq = J8aG Tpye ii, C2) ceay af (A2.05 ) 


0 


I eA 
Zan 2 agg | sin® (PL) L) Ot [x2 5? =) df 
a ¥ (Py)? free?’ 1, 62)( fac? < a) 

O 


7 | 
J sin2(PL) (+3) Ko( VP 2k? a) a 
Con (-g) Pek? Ky ( [7 2-2 a) 


(A207) 


t 80% 
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To reduce (A2.07) to the form of Robey's equation (18) consider the 


following integrals: 


‘i ye 
ean () L) 


Coe 


Ii = 


Ieee J = ke sini 

then when J = 0, &= 0 
yek, ws 7 

end \k2 -/2 - oe 


dfeaekcoswdw 


Therefore 
ar 
2 
is aan Cee sint? ) 
tn? aS 
O 5 
Consider 
00 
_ | sin*(P L) et 
i= = 
(J) 
k 
Let P= k coshy 
then when f~ =k, #2 0 
%= 00, V= 


! 
and | pak? =k sinhY 


df = k sinh dy 
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He?) (ke cos 7) 


Hy \?) (ca cos wt) 


1 
eel a, ) 


2.46 Ky ( penne! a) 


» 


(AZ 00m) 
a= (A2.08 ) 
a? (A2.09) 





Therefore 


00 
ae) i 
ee a: sin“(kL cosh¥ ) Ko (ka sinh¥Y ) ay (A2,10) 
Ke (cosh ¥ Me ; 
0 Ky (ka sinh ¥ ) 
Equation (A2.05) may now be written as 
a 
\" ( 
aq = 5 8a PoW | sin = aL, sw) 2) (ice cos 7) que ws 
~ Z 
. 0 (sinw)” rt? (ee cos & ) 
1, 
re sin” (kL se cosy K 9 (ka sinh} ) | 
A oan e K, (ka sinh / ) 


Thus from (A2.11), it can be seen that the terms in Robey's equation 
(18) should have identical Signs; and from the initial argument about 
equation (16),- the terms should be positive. 

Robey's expressions for fr end as they appear in equations (19) and 
(20) will now be derived from (A2.11). The C torm is self explanatory and 


correct. Examine the expression 


a 2 
Hp (2) (ka cos  ) = Ho (2) (x where x = ka cos w 


H, $2) (ca cos tw ) H, (?? (x ) 


Jo(x) = j¥o(x) 
Jz (x) a JY (x) 
Jo({x) = j¥aG Ty (x) Yz (x) 


3 (x) Zvi (x 
- Jo(x)Jy (x) #Y,(x)Yz (x) Bees 1 (x) ~ Jy (xe ee 


Fix) AYE (x) Jy (x) A ¥y (x 
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¢ 


Jo(ka cos 2*)J ,(ka cos 2») AYo(ka cost) Yj (ka cos 2) 


z 2 
J, (ke cost) FY (ke cos t") 


e 7 
io ee l 
ka cos wv 





Za 
J, (ka cos @%) AX (xe cos z~) 


Thus it can be seen that Robey's coefficients “&,43 , and € in 
equations (19) and (20) are correct as stated but are not consistent with 
equation (18). 

The above arguments are believed to be sufficient to show Robey's 
equation (18) is not correct as stated, but would be correct in the form 
of equation (A211), Equation (A2.11), if stated in terms of average 
radiation impedance and with a change of variable, is equivalent to 


equation (3.02). 


to, N. Watson, A Treatise on the Theory of Bessel Functions, The Macliillan 
Company, 1945, p 77 
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